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LONG-TERM  GOALS 

Random  variability  in  shallow  water  will  induce  variability  in  a  propagating  acoustic  field.  The  long¬ 
term  goal  of  this  research  is  to  quantify  how  random  variability  in  the  ocean  environment  translates 
into  random  variability  in  the  acoustic  field  and  the  associated  signal  processing  algorithms  in  the  mid¬ 
frequency  (1-10  kHz)  band.  In  the  present  funding  cycle,  the  emphasis  is  on  the  waveguide  invariant. 

OBJECTIVES 

Constructive  and  destruction  interference  is  an  inevitable  consequence  of  multipath  acoustic 
propagation  in  shallow  water.  If  the  ocean  environment  is  sufficiently  benign,  the  so-called  waveguide 
invariant  describes  the  resulting  interference  pattern.  The  waveguide  invariant,  usually  designated  as 
“beta,”  has  traditionally  been  regarded  as  a  low-frequency  phenomenon.  The  objectives  of  the  present 
work  are  twofold:  first,  to  extend  the  waveguide  invariant  concept  so  it  may  be  incorporated  into  mid¬ 
frequency  signal  processing  algorithms,  and  second  to  quantify  the  limitations  on  these  algorithms 
imposed  by  shallow-water  internal  waves. 

APPROACH 

The  approach  is  a  mixture  of  data  analysis,  theoretical  development,  and  numerical  modeling.  Dr. 
Dajun  Tang  of  the  University  of  Washington  Applied  Physics  Laboratory  (APL-UW)  is  a  key 
individual  in  acquiring  suitable  acoustical  and  environmental  data  sets  that  can  be  used  in  the  analysis. 
Dr.  Lisa  Zurk  of  Portland  State  University  is  a  key  individual  in  developing  the  signal  processing 
algorithms.  Dr.  Andrey  Lunkov  of  the  General  Physics  Institute,  Moscow,  has  been  a  key  individual  in 
developing  theoretical  models  for  calculating  array  coherence.  Dr.  Lunkov  received  support  for  the 
collaboration  from  ONR-Global. 

WORK  COMPLETED 

Under  previous  support,  the  author  developed  a  signal  processing  method  for  calculating  the 
waveguide  invariant  distribution  [Rouseff  and  Spindel,  2002],  In  this  approach,  rather  than  try  to 
assign  a  single  value  to  the  waveguide  invariant,  a  distribution  of  possible  values  is  calculated.  The 
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distribution  may  be  broad  or  sharply  peaked  depending  on  the  environment.  Numerical  simulations 
showed  how  shallow  water  internal  waves  could  broaden  the  distribution  and  possibly  shift  the  peak 
[Rouseff,  2001].  The  waveguide  invariant  was  incorporated  into  signal  processing  algorithms  for 
coherent  processing  of  array  data.  A  method  called  striation-based  beamforming  was  developed  to 
improve  array  gain  [Rouseff  and  Zurk,  2011;  Zurk  and  Rouseff,  2012],  It  was  also  shown  theoretically 
how  striation-based  beamforming  permits  one  to  estimate  both  the  source-to-receiver  range  and  the 
waveguide  invariant;  it  had  previously  been  thought  that  only  the  ratio  between  the  two  quantities  was 
accessible 

Under  current  support,  striation-based  beamforming  was  tested  using  data  collected  during  the  201 1 
Gulf  Experiment  (GulfExl  1).  Both  the  source-to-receiver  range  and  the  numerical  value  of  the 
waveguide  invariant  were  estimated  every  0.5  Hz  over  a  200  bandwidth  centered  at  3.5  kHz.  Results 
were  presented  at  the  2nd  International  Workshop  on  Acoustic  Interference  Phenomena  and  Signal 
Processing  [Rouseff,  Zurk  and  Tang,  2013],  The  workshop  was  held  in  Moscow  in  2013  with  the  PI 
serving  as  co-chair  and  co-organizer. 

Efforts  in  FY14  returned  to  studying  how  shallow-water  internal  waves  affect  the  waveguide  invariant. 
An  analytical  expression  was  derived  for  the  average  waveguide  distribution  in  the  presence  of  random 
internal  waves  described  by  a  modified  Garrett-Munk  spectrum.  Analytical  expressions  were  also 
derived  for  he  expected  array  coherence  length  using  the  same  statistical  model  for  the  internal  waves. 

RESULTS 

Figure  1  is  a  sample  calculation  showing  how  shallow  water  internal  waves  cause  coupling  between 
acoustic  modes.  The  plot  shows  the  intensity  of  the  propagating  acoustic  modes  as  a  function  of  range 
for  a  typical  shallow  water  scenario.  The  starting  field  is  the  eighth  acoustic  mode,  a  mode  that  is 
bound  from  above  by  the  sea  surface.  As  the  range  increases,  energy  leaks  out  of  the  eighth  mode  into 
lower  order  modes.  These  lower-order  modes  are  bound  from  above  by  the  thermocline  rather  than  the 
sea  surface.  The  result  is  significant  as  refracted  lower-order  modes  have  a  different  value  for  the 
waveguide  invariant  than  higher-order  reflected  modes.  Consequently,  at  sufficient  range,  the  apparent 
value  for  the  waveguide  invariant  will  change. 

Figure  2  is  a  sample  calculation  showing  how  shallow  water  internal  waves  can  change  the  waveguide 
invariant  distribution.  The  calculation  uses  the  analytical  expression  derived  during  the  current  fiscal 
year.  Internal  waves  broaden  the  distribution  and  shift  the  peak  to  larger  values  of  beta. 

The  calculations  in  Figures  1  and  2  were  performed  assuming  acoustic  frequency  400  Hz.  At  present, 
we  are  working  to  extend  the  calculations  to  the  mid-frequency  (1-10  kHz)  regime.  These  extensions 
are  necessary  as  the  acoustic  waveguide  will  support  many  more  modes  at  the  higher  frequency  and, 
without  modification,  the  calculation  becomes  cumbersome.  The  analytical  predictions  can  then  be 
tested  using  towed-source  data  collected  during  the  2006  Shallow  Water  (SW06)  experiment. 

Results  derived  during  FY14  have  been  presented  at  two  conferences  [Rouseff,  2013;  Funkov  and 
Rouseff,  2014],  Efforts  in  FY15  will  include  documenting  these  results  in  refereed  journal 
publications. 


2 


1 


.§■  01 

m 

c 

s 

e 

QJ 

"Cl 

Q 

E 


range  (km) 


-13 


-10 


Figure  1:  Average  modal  intensity  as  a  function  of  range .  Plot  shows  how  energy  that  starts  in  mode 
8  leaks  into  other  modes  as  the  range  increases  due  to  ocean  internal  waves.  The  calculation  is  for  a 
typical  shallow-water  environment  in  the  summer.  Results  are  for  400  Hz. 
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Figure  2:  Average  waveguide  invariant  distribution.  Calculation  for  a  source  in  the  surface  mixed 
layer  and  a  receiver  20  km  downrange  and  below  the  thermocline.  Plot  shows  internal  waves 
broaden  the  distribution  and  shift  it  to  higher  values  of  the  waveguide  invariant,  beta. 
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IMPACT/APPLICATIONS 


The  waveguide  invariant  is  currently  applied  in  practical  sonar  signal  processing  algorithms.  Extending 
its  usages  to  the  mid-frequency  regime  is  desirable. 

RELATED  PROJECTS 

This  project  uses  acoustical  and  environmental  data  collected  in  ONR-supported  experiments  like 
SW06,  GulfEx  and  TREx.  Collaboration  with  investigators  from  APL-UW  in  using  these  data  sets  will 
continue.  Collaboration  with  investigators  from  Portland  State  University  supported  by  the  ONR 
Underwater  Signal  Processing  will  also  continue. 
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